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As a first step toward the understanding of Rayleigh-Brillouin scattering from a fluid with a temperature 
gradient, we analyze the initial value problem for certain prototype one-dimensional nonuniform systems. For 
sufficiently short times and localized initial pulses, it is not necessary to impose actual physical boundaries on 
the linearly nonuniform system. By a straightforward, though unconventional, application of Fourier and 
Laplace transform methods, explicit and physically reasonable solutions are constructed for the propagation of 
fluctuation pulses in space and time according to the nonuniform wave and diffusion equations and a special 
case of the nonuniform damped wave equation. The analog of the dynamic structure factor is also constructed 
for the latter two cases. 

Key words: damped wave motion; density fluctuations; diffusion; dynamic structure factor; initial value problem; 
nonuniform systems; one-dimensonal models; Rayleigh-Brillouin scattering. 

I. Introduction 

There has been considerable recent interest in the problem of Rayleigh-Brillouin scattering in a fluid with 
a stationary temperature gradient. [1-10] 1 The problem, in essence, is to generalize Mountain's standard 
derivation of the equilibrium Rayleigh-Brillouin spectrum [11] to the nonequilibrium stationary state. 

Laser light is scattered from a liquid because of random density (and therefore dielectric constant) fluc- 
tuations. [12, 13] In the hydrodynamic regime, a density fluctuation will propagate in time partly as a diffu- 
sion or "heat" mode, causing the Rayleigh line, and partly as a wave or "sound" mode, causing the two 
Brillouin lines. With a temperature gradient the spectrum is altered. For example, if more sound modes are 
generated at the hot side of the liquid than at the cold side, there will be a larger flux of sound modes from 
hot-to-cold than from cold-to-hot, thereby causing a size asymmetry in the Brillouin peaks. Some treatments 
of this problem, [1-4] in fact, predict no spectral changes other than an overall asymmetric scale factor in 
the Brillouin peaks. Such an asymmetry has recently been observed experimentally. [14] 

Other authors [5-9] predict changes in the mathematical shape of the Brillouin lines from their 
equilibrium (Lorentzian) form. Such changes of shape, very difficult to measure experimentally, are inter- 
preted here as the effect of propagation of modes through regions of changing sound speed, sound attenua- 
tion, and thermal diffusivity. Among the theories published to date, [1-10] there is significant disagreement 
as to the precise form of the nonequilibrium spectrum. Some discrepancies may arise from different assump- 
tions about the nature of the nonequilibrium steady state, but it is not clear at this time if all do. 

None of the theories published thus far appear to follow the procedure of ref. 1 1 explicitly, i.e., to solve 
the hydrodynamic equations for the time evolution of a density fluctuation and, subsequently, to perform a 
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stochastic and ensemble average to determine the dynamic structure factor and hence the spectrum of scat- 
tered light. Closest in spirit to such a procedure is the hydrodynamic approach of van der Zwan and Mazur, 
[7, 8] but they appear to bypass the calculation of the time development of an individual fluctuation. 

An explicit calculation of the evolution in time and (both physical and Fourier) space of a spontaneous 
density fluctuation, in the presence of a temperature gradient, should shed considerable light on the 
presently confused predictions of Rayleigh-Brillouin scattering. The possibility of changes with tempera- 
ture, and hence distance, of all physically relevant variables should be considered. Significantly, possible 
changes in the shape of the Rayleigh line have been substantially neglected thus far. [1-10] 

This paper represents an initial step in the execution of the above program. Since the full three- 
dimensional problem of coupled equations in density and temperature fluctuations is very complex 
mathematically, it is illuminating to deal first with simplified models. Here we solve explicitly the one- 
dimensional wave, diffusion, and damped wave equations to first order in nonuniformity. Care is taken to 
select nonuniform equations which correspond as closely as possible to the physical nonequilibrium steady- 
state problem. 

Our technique is a generalization of a method due to Brownell [15] to calculate the Green's function for a 
nonuniform wave equation. We show that the initial value problem for the nonuniform wave, diffusion, and 
damped wave equations can indeed be solved, to linear order in nonuniformity, in a manner which closely 
parallels Mountain's derivation. The methods are systematic, but the mathematics becomes rapidly more 
difficult with greater complexity of the underlying differential equations. Our eventual goal, yet to be 
attained, is the explicit solution of the coupled hydrodynamic equations and determination of the dynamic 
structure factor in one and, if possible, three dimensions. Since wave and diffusion equations are prevalent 
throughout physics, the techniques presented here may prove useful for a variety of other physical problems 
as well. 

We begin in section 2 with a general discussion of some predictable mathematical difficulties. Section 2.1 
examines the applicability of the initial value problem to Rayleigh-Brillouin scattering in a fluid with a 
temperature gradient. Section 2.2 describes the needed transformations of variable and examines the range 
of validity of a first-order nonuniformity calculation. Solutions to the various initial value problems in the 
uniform limit, mostly well-known, are outlined in Section 3. The subsequent three sections present explicit 
solutions, with intuitively oriented interpretations, to the initial value problem for the nonuniform wave, dif- 
fusion, and damped wave equations respectively. Results are summarized in Section 7. 



2. Initial mathematical considerations 
2. 1 . The initial value problem 

The basic relationship between the measured spectrum of scattered light and the structure of a liquid 
sample is [12, 13] ^ 

/(w) = C S (k, w) (1) 

where co is the shift in frequency and k the shift in wave vector of the scattered light, / is the average 
measured intensity , and C is a constant, not important for our present interests. 5(k, <u), the dynamic struc- 
ture factor, is the Fourier transform in time of the intermediate scattering function F(k, 

5{k,w) = ]Z e-«" F(k, t) dt (2 ) 

where 

F(k,<)= <H(-k,0)«(k,0> - (3) 
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Here the brackets denote an ensemble average and n (k, t) is the Fourier transform in space of n (r, t), the 
fluctuation of liquid density about equilibrium, i.e., 

n (k, = j dr e- ,kr n (r, t) . (4) 

The hydrodynamic regime is defined by the inequality 

k- 1 > > I (5) 

where I is the interparticle spacing, and 

*=2*,sin^-, (6) 

where k, is the wave vector of the incident light and $ is the angle by which the light is scattered. In a liquid t 
is small and scattering experiments with visible light are always in the hydrodynamic regime. For dilute 
gases ( is large, but the hydrodynamic regime can be attained [16] with sufficiently small 6. Kirkpatrick, et 
al. [5] consider the nonuniform problem starting from the Boltzmann equation, which is strictly valid only 
for dilute gases. The initial experimental observation of changes in the scattered spectrum due to a 
temperature gradient has been made with liquid water at a very low scattering angle. [14]. 

Mountain's derivation [11] of Rayleigh-Brillouin scattering in equilibrium begins with the observation of 
Landau and Placzek [17] that, in the hydrodynamic regime, density fluctuations, n(r, t), and temperature 
fluctuations, 7\r, f)» satisfy a set of coupled second-order linear differential equations. After a Fourier 
transformation in space and a Laplace transformation in time, those equations become coupled algebraic 
equations which involve the initial {t = 0) values and time derivatives of the fluctuations. The derivation of 
ref. [1 1 ], in effect, makes the further simplifications 

n (k, 0) = (7) 

f(k, 0) = (8) 

where a dot denotes a time derivative. These simplifications are not literally correct, but are permissible 
because, at the end of the calculation when the autocorrelation function is evaluated, the above quantities 
are not correlated with the initial density fluctuations. In particular, because of time reversal invariance, 

<n(- k,0)n(k,0)> =0 (9) 

and because of the thermodynamic independence of density and temperature, 

<n(- k,0) r(k,0)> = 0. (10) 

The solution of the differential equations leads to a dynamic structure factor, eq. (2), in terms of the static 
structure factor 

S(k) - <n(- k,0)n(k,0)> (11) 

which must be determined from separate (but well-known) thermodynamic considerations. [12, 18] An alter- 
nate method to the above initial value problem is to consider the response of the system to a random fluc- 
tuating force. [19] 

In the presence of a stationary temperature gradient, time-reversal symmetry no longer holds, [1, 4] and 
eq (9) is no longer satisfied. Thus an initial value calculation, by itself, does not suffice to determine the 
altered Rayleigh-Brillouin spectrum. However, an initial value approach may suffice to determine "shape 
effects," as opposed to "size effects," in the altered spectrum. 
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More precisely, we define e, which is proportional to | VT | , as the appropriate dimensionless expansion 
parameter. The LHS of eq (9) will be 0(e) due to differences in the rate of generation of fluctuations at dif- 
ferent positions in the fluid. These differences will lead to overall asymmetries in the Brillouin lines. A quite 
different effect is the change in the shape of the Rayleigh and Brillouin lines due to propagation of the 
unperturbed fluctuations through a nonuniform medium. 

In the linearization of the problem about uniformity, it is reasonable to assume that the two effects are 
independent and additive. Further corrections due to the nonuniform propagation of fluctuations which are 
themselves generated by nonuniformity presumably lead to second order, 0(€ 2 ) corrections. 

Thus an initial value problem approach to Rayleigh-Brillouin scattering in a temperature gradient, while 
incomplete, may suffice to determine the shape changes of the lines, and hence will complement those 
approaches which appear to emphasize only size changes. [1, 4] 

2.2. Transformations and nonuniformity expansions 

In the presence of a temperature gradient, any relevant physical property A of a fluid may be expanded in 
a Taylor series about some convenient origin of position (r = 0) in the fluid, i.e., 

A{v) = ,4(0) + (r • VT) {-jfirU + . . . . (12) 

For the linearized problem, only the above two terms are retained. 

The coefficients of the hydrodynamic equations, in equilibrium eqs. (6)-(8) of ref. [11], will then be 
position-dependent according to eq (12). Since the differential operators do not commute with r, it is impor- 
tant to return to the fundamental derivations of the differential equations governing the processes, to deter- 
mine in what order differential operators and position-dependent coefficients should be placed. 

In the prototype one-dimensional equations considered here, the coefficients will have the position 
dependence 

a(x) = a + a\X (13) 

where d t is assumed to be small. 

We will examine for the various differential equations the initial value problem, i.e., the solution for y(x,t) 
given y(x,0) and zero initial time derivative, 

y (x,0) = . (14) 

Equation (14) is not necessary for the diffusion equation, which is first order in time. 

It is convenient at the outset to introduce a special notation for the several variable spaces which are re- 
quired. Physical space or "P-space" denotes the variables [x,t\. Fourier space or "F-space" consists of the 
variables [k,t], where 

y(k,t)= n. dx e-<* y(x, t) . (15) 

Laplace space or "L-space" is formed by the variables [x,s], where 

y(x,s)= L". die'" y(x, t) . (16) 

Finally, Fourier-Laplace space or "FL-space" consists of the variables [k,s], where 

y\k, s) = il dx e -«- y(x, s) = L". dt e- y(k, t) . (17) 

Our convention for transformed functions is that a bar denotes F-space, a caret L-space, and a tilde 
FL-space. 
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The differential equations considered here require that certain coefficients be positive, which is, in 
general, in conflict with eq (13). For small but finite a„ if 

x < — aja.1 

then a(x) < 0, which could lead to physical nonsense, e.g., negative diffusion or imaginary wave speed. We 
must assume that such large negative values of position are irrelevant to the problem. 

A possible escape from this difficulty is to impose explicit boundaries on the linearly nonuniform region 
of space with, for example, separate uniform regions outside those boundaries. 

a(x) = a (,) x < - I 

= — [a (1 > + a' 2 '] 
2 J 

+ _L[ a u> _ a u>]JL _ L < x < i 

= a< 2 > x > I 

We elect, however, not to proceed in this manner because a new parameter L is introduced into the prob- 
lem, and the exact solution for a general L is likely to be mathematically cumbersome. Furthermore, for the 
conditions we seek to impose, there should exist a unique linearized solution 

y(x, t) = y (0) (x, t) + a,y (,) (x, t) + (Xe 2 ) , (18) 

where y <0) is the uniform solution, which is independent of L. These conditions are that the initial value 
packet be localized and that the time be small compared to the time required for the packet to propagate (in 
significant measure) to the region of negative a. We consider initial value packets of finite width w , i.e., 

y(x,Q) = \ x \ >f (19) 

and the "smallness" of aj is made precise in dimensionless form by 

2iBl < < l . (20) 

Go 

Additionally, we require a constraint on the time t of propagation. For example, the solution to the 
uniform wave equation moves with constant velocity c, so we require times short enough that 



aid 



<< 1 (21) 



or, for the diffusion equation, since the mean-square width of an initial delta-function pulse is (in one 
dimension) 

<x 2 > = 2Dt (22) 

the appropriate condition on the time is 

a, (2Dt) u2 /a < < 1 . (23) 

Calculation of the dynamic structure factor, in principle, requires an integral from t = to t = oo. Never- 
theless, we expect the linear nonuniformity theory to give a meaningful structure factor if the fluctuation 
has substantially decayed away in a time such that the inequalities (21) and (23) are valid. The phrase 
"decay away" is imprecise; a norm of the fluctuation needs to be defined and must be shown to be negligi- 
ble after some particular time interval. 
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One possible norm is 



\y\ - Too dxy(x, t), 



(24) 



but, for the equations considered here, eq (24) happens to be a constant of the motion, and does not decay in 
time at all. 

What do decay in time however, are Fourier components of the fluctuation for finite h. The h = com- 
ponenet is \y \ itself. For diffusion, y\k, t) decays as e~ Dl ^', and for the damped wave equation considered 
here, it decays as e" r * ', where T is the attenuation coefficient. In the Rayleigh-Brillouin experiment, a par- 
ticular value of h is monitored by virtue of the fixed laser light frequency and scattering angle. If we are 
similarly interested in monitoring a particular value k = k , then for the times 



t D = (d kiy 



(25) 



and 



ip = (r kir 



(26) 



we require that inequalities (21) and (23) hold. The various inequalities for our particular applications in 
terms of variables defined in later sections are listed in table I. 



Differential 
Equation 



Table 1. 



Conditions for 
Validity of 



Conditions for 

Validity of 

S(Jfc,u>) 



Undamped wave 



bw n 



bet 



<< 1 



<< 1 



Never valid 



Diffusion 



bw n 



<< 1 



ak 



<< 1 



Damped wave 



-(ZV) ,/2 << i 

a 
"- << 1 



rt 



be 



<< i 



(_LL <<1 ) _^£L 



<< l 



The standard method of solution to the initial value problem for a partial differential equation in P-space 
is to transform to FL-space. In the uniform cases, the equation is an algebraic one and is easily solved. The 
solution is completed upon inverse transformation to P-space. In ref. [11] only the inverse transformation to 
F-space is required to obtain the intermediate scattering function, eq (3), and the dynamic structure factor, 
eq(2). 
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Since, under Fourier transformation 



dk 



a partial differential equation in P-space with linearly nonuniform coefficients becomes, in FL-space, a 
linear differential equation in k. It is always possible to solve a linear differential equation by the integrat- 
ing factor technique. [20] However, a unique solution requires a boundary value on k, and the problem 
appears at First not to specify any particular boundary value. We shall see that this apparent new degree of 
freedom is bogus and that the physics imposes a particular boundary value on k. 

A more cumbersome problem is the inverse transformation. As stated previously, in P-space the solution is 
nonzero only within a width w about the origin such that (20) holds. In F-space, according to Heisenberg's 
inequality, [21] the solution y{k, t) has a width u£'. Only the values of x such that 



< < 1 



(27) 



are significant; however, in F-space there is no analogous inequality. For example, 

<< 1 



a e \k\ 



(28) 



does not hold for all significant values of \k\, since such values of ] k \ can range from zero to w~ B v . 

We conclude that there is no natural means, in FL-space or F-space, to take advantage of the smallness of 
the nonuniformity parameter. Therefore, after solving the linear differential equation in FL-space, we must 
first inverse transform to i-space (and expand in nonuniformity where convenient), then inverse transform 
to F-space (and further expand in nonuniformity), and lastly transform to F-space. This circuitous but 
required route is shown schematically in figure 1. 

Expand in nonuniformity 



P space 

y(x,t) 

I 

i 

I 
F space 

y(M) 

I 

i 
I 

\ 
F(k,t) 



L space 

y(x,s) 

♦ 

i 

i 

FL space 

y(k,s) 



S(k,u) 



FIGURE 1. Sequence of transformations for the uniform problem (solid ar- 
rows) and for the nonuniform problem (dashed arrows). 

In the following sections we carry out the above procedure explicitly for three nonuniform differential 
equations, after first reviewing their respective uniform solutions. 
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3. Uniform differential equations 

In this section, for reference purposes, we briefly outline the soluions to the initial value problem for the 
one-dimensional uniform wave, diffusion, and damped wave equations. These results are well-known. 

3,1, Wave equation 

The wave equation in P-space is 

_^y__J_J!2L = , (29) 

dx 2 c 2 dt 2 

We assume a given initial wave packet y{x, 0) with zero initial time derivative. The equation is FL-space is 

(30) 



s -/ 



- k 2 ftk, s) - i- y\k, s) + -^y(k, 0) = 
c c 



with the solution 



The solution in F-space is 



the solution in Z-space is 



= Vfo 0) (31) 



R k ' s) " s 2 + c 2 A 2 * 



y(*» = 4- yU. °) [e fcte + e ~ tM \ > (32) 



f{x, s ) = _L £ da e-'\"\ ,c y{x - u,0) (33) 

2c 



and the solution in P-space is d'Alembert's solution, 

y(x, t ) - —■ \y{x - ct, 0) + y{x + ct, 0)], (34) 



2 



i.e., the pulse divides in half and each half moves, without change of shape, at speed c. For an ensemble of 
packets, the analog of the intermediate scattering function is 

<y{- k, 0) y\k, t)> = 4" <y(~ K 0) y{k, 0)> [e Ukt + e' Uki \ (35) 

2 

and its Fourier transform in time, the analog of the dynamic structure factor, is 

S(k, u) = it <y{-k, 0) y(k, 0)> [fi(u - ck) + ofa + ck)], (36) 

i.e., two "Brillouin lines" of infinitely small width, since there is no damping. 

3.2. Diffusion equation 

The diffusion equation in P-space is 
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An initial wave packet y(x, 0) is given. Since the equation is first-order in time, no initial time derivative is 
specified. The equation in FL-space is 

syik s) - y(k, 0) + Dk 1 y\k. s) = (38) 

with solution 

The solution in F-space is 

y\k, t) = y{k, 0) e- Dkh , (40) 

the solution in L-space is 

y(x, s) =±.( s Dr» H dvy{v, 0) e^^'V-l, (41) 

and the solution in F-space is 

y(x, t) = 2( T j f) «/i 11 dv y(v, 0) e-i* - •>*'«» . (42) 

For a 5-function initial condition, y{x, 0) = 8{x), the solution takes the familiar form 

**-» = ^rw^ <'■*""■ (43) 

The analog of the intermediate scattering function is 

<y(- k, 0) y{k, t)> = < y{ - k, 0) y(k, 0) > e" D * 2 « (44) 

and the analog of the dynamic structure factor is 

S(k, a) = < y{- k, 0)y(k 0) > ^ f^frf • < 45) 

i.e., a Lorentzian "Rayleigh line" with width Dk 2 . 

3,3. Damped wave equation 

There are several possible choices for a wave equation with damping. We choose the case that cor- 
responds most closely to the Rayleigh-Brillouin problem, i.e., an equation with "P-damping." In P-space 
this is 

JV + JLJi.JL_lA = (46) 

dx* c 2 6V dt e be u " yw) 

We again begin with an initial wave packet y{x, 0) with zero initial time derivative. The equation in 
F-space is 

k*y(k, t) + ^-^yik, t) + -1 ^y(k, t) = (47) 

which is the equation for a damped harmonic oscillator. The condition that the oscillation be underdamped, 
which we impose here, is 

Vk 

-^- < 1 . (48) 

c 

437 



The equation in FJL-space is 

A 2 ftk. s) + ^- y[k, s) - -^/- y(k, 0) + £ y\k, 5) 



--£y(A.0) = , (49) 



with the solution 

s + 2Tk 2 
y[k, s) = y{k, 0) 7-^2nJFT^F 



(50) 



The solution in F-space is 

y(k, t) = -i-flfcO) {e** [e>"*-r»>>">» + o-.^-rV,"^] 

+ ^-rW" *""'' [*"• 1 -'• 2 ■ ,, " 2 " " .-^"''"'l} (51) 

and the solution in L-space is 

** -) = ^rkp- '- dv "** 0) ~ "T" ^ M)1 

Xe-l~l'». a W a . (52) 

However, to the best of our knowledge, there is no representation of the solution in P-space in terms of 
elementary functions. But a convenient approximate solution in P-space may be obtained by means of a 
dispersion expansion. First, the concept of a A-dependent wave speed is introduced, i.e., 

e(A) = c (l-I^l)-. (53) 

If -T± < < 1 (extreme underdamping) we may expand in this small "dispersion parameter." Then eq (44), 
to leading order, is 

y(k, t) = ~y(k, 0) {e-" 2 ' [<?'«*' + e""*']} (54) 

with a corresponding solution in P-space 

yix, t) = ^ t)l/2 II dvy(v, 0) [ e -<-— > a '«* + e -^«^<] , (55) 

which corresponds to two wave pulses travelling in opposite directions and spreading with time in a diffu- 
sionlike manner. In the nonuniform case there will be expansions in both dispersion and nonuniformity; 
however it is reasonable to assume that the two expansions are independent of each other. 
From the lowest-order solution, eq (54), the analog of the intermediate scattering function is 

<y\- k, 0) y(k, t)> = ~ <y(- k, 0) y(k, 0)> e'" 2 ' [e ukt + e~ ickt ] (56) 

and the analog of the dynamic structure factor is 

* -» = <*- fc 0) **• 0)> x K-JHwvt + ^IyHV^ • (57) 

i.e., two Lorentzian "Brillouin lines" each with width Yk 2 . 
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4. Nonuniform wave equation 

Having reviewed the initial value problem for the three uniform differential equations, we next derive 
solutions to their nonuniform counterparts. The solutions must have the form of the uniform solution plus a 
small correction linear in nonuniformity. In all cases the nonuniformity parameter will be denoted by the 
variable b. The solutions are subject to the constraints on time, distance, and wavenumber described in sec- 
tion 2.2. 

In the following three sections we first present the exact form of the nonuniform differential equation 
which we take to be (in one dimension) closely analogous to the three-dimensional Rayleigh-Brillouin experi- 
mental situation described earlier. Thereafter we present a detailed mathematical derivation of the solution. 
Each selection then ends with a physical interpretation of the solution thus derived. 

4.1. The differential equation 

Following Brownell, [15] we write the nonuniform wave equation as 

Q (*)-^-T^- = (58) 

where q(x) = a + bx. (59) 

Physically, this corresponds to a stretched string with a variable mass per unit length q(x) at a constant 
tension T. This is analogous to the fluid with a stationary temperature gradient, where the pressure is con- 
stant throughout but the density varies (to first order linearly) with distance due to the gradient in 
temperature. Again we consider a given initial pulse y{x, 0) with zero initial time derivative. 

4.2. Solution 

The equation in FL-space is 

The first-order differential equation in k is solved by standard methods [20] as 
y[k, s) = \i dk'[-£-M 0) - j-lrjli; 0)] 

Xexpl - b ^ J. (61) 

Formally, the integration constant k x appears to give a new and unwanted degree of freedom to the prob- 
lem. However, as pointed out by Brownell [15], the existence of an inverse transformation requires certain 
analyticity properties in the complex s-plane. In our formulation, the requirement is that y{k, s) be 
analytic in the right half of the complex s-plane, which in turn requires that (k' 3 -k 3 ) be the same sign for all 
it. This requires that it, = oo. We also henceforth assume, without loss of generality, that b is positive. 

Changing variables to u = k'-k yields 

y[k, s) = %du [-£-*« + a. o) - -j-Hryi" + k > °J 

x ex P [~ "ir - 357 "(" 2 + Uu + 3 ^] ' (62) 
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As explained in section 2.2, we must now transform to £-space. Let 9" 1 denote the inverse Fourier 
transform operator, i.e., 

9" l/W) = 4r I: « dk etk *M • (63) 

Now, 

9- \-^-y(u + k, 0) - — -|-y(« + h 0)1 

= ^-y(^0) e — (1 +^). (64) 

os a 

Next we require the inverse Fourier transform of the second term in eq (62). The relation 

j- dk e -^-Bk = (JL) ,/2 e e 2 «« (65) 

holds for complex a and providing Rea > 0. To apply eq (65), we must assign a small negative imaginary 
component to T, i.e., 

T-T- ie (66) 

with the limit e — taken at the end of the calculation. It then follows that 

9-1 {exp [- Jf- - JL u(u> + Mu + 3i»)]} 

= i- &huiT?» exp [± (-££- - 2«* + ^)] (67) 

[ma Tiu 3 1 

" 1 3fcH ■ 

The assignment of a small imaginary part to the coefficient T appears arbitrary. It may be justified 
mathematically by an extension of BrownelPs method, [15] which uses the complex Laplace transform with 
real coeffients and appropriate contours off the real axis in the complex s-plane. Here we employ the opera- 
tionally simpler, and presumably equivalent, method of complex coefficients and a real Laplace transform 
variable. 

From eqs (62), (64), and (67), and the convolution theorem for Fourier transforms, we obtain 

y{x, s) = — ai (iriT)~ in \™dw w~ U2 exp (- law — — ^ ) 

X J". dvy{v,0)exp(-ibvw){l + — ) (68) 

a 



, af {i[i^-^-„ + « 



Tw 

where w = bv. At this point we may discard terms which are explicity 0(6 Z ) and rearrange, with the result 
that 

Six, s) =~ai(iviT)- ul I.". dvy{v, 0X1 + — ) 
z a 

X %dw w- in exp[-iaw - -~ bu{x + v) + is 2 (x - vfllTw], (69) 
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To evaluate the w-integral, we make the substitution w = z 2 . The relation [22] 

g-e-W dz = ^ (jl)"> e -2c^,«« (70) 

holds for all a, such that /?e(a) > and i?e(/3) > 0. To apply this relation to eq (69), in addition to the 
small negative imaginary term added to T, a small positive imaginary term must be added to a. Then eq (69) 
reduces to 

y{x, s) = \{a/T)*>* £. dv y(v, 0)(1 + — ) [1 + £(x + t;)]"" 2 

X exp \- $ \x - v\ [(a + yfe [x + V J)IT] U2 } . (71) 

At this point it is evident, by inspection, that the uniform limit b — + recovers the uniform solution, eq 
(33), where c = (T/a) in . 
To find the P-space solution, we use the relation 

£-\ | e -o,j = 5( t - a ) (72) 

where-£ is the Laplace transform operator. This holds for real positive a, and hence is applicable provided 
that 

(* + v)> -y. (73) 

Thus we see that the concerns raised in section 2.2 about the absence of explicit boundaries are built into 
the mathematics of the solution; an inverse Laplace transform cannot be performed for large negative x, 
where the differential equation becomes physical nonsense. 
With this proviso, the P«space solution is 

yix, t) = \{a/n n L". dv y(ifr0)(l + y) [1 + A (* + „)]-"» 

%t-\x-v\[a+^b{x + v)] tf2 IT) (74) 

which, to leading order in nonuniformity, may be reduced to the more suggestive form, 
y(x, t) = -1(1 - -—-)}<* ~ ct [1 - A (2x - ctm 

+ jo + -%*-)& + c *p -h®* + c0] » 0) (75) 

where c = (77a) 1 ' 2 , i.e., two moving pulses. It is easy to show that, through terms linear in b, the pulses 
separately satisfy the original differential equation, eqs (58) and (59). 

The F-space solution may be derived from eq (75) and, to leading order in nonuniformity, can be 
represented in a similarly suggestive form, 

y{k, t) = i-(l - ^-)y{k[\ - -^-l,0)«p[-ic(l - ~^ht\ 

+ T (1+ 1T ) ^ [1+ 4f ] ' 0) exp[ ic (I + "^ kt] ' (76) 

cf. eq (32). 

Since eq (76) holds only for small times (cf. table 1), and there is no damping for the unperturbed solution, 
it is not meaningful to construct a dynamic structure factor for this case. 
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4.3. Interpretation 

It is at this point easily seen that the solution has the proper physical behavior expected on intuitive 
grounds. The first term of eq (75) is the rightward moving pulse. Since b is positive, the wave speed 
decreases to the right, and, as it should, the pulse moves more slowly as time progresses. The factor 

(1 ££L) indicates that the pulse size decreases with time, as it should since it moves into a region of 

4a , 

greater inertia. In addition, the width of the pulse gets smaller in time by a factor (1 + — — )~\ because the 

front end of the pulse moves more slowly than the back end. Conversely, the leftward moving pulse, or sec- 
ond term of eq (75), gets larger and wider and moves faster with time. 

The F-space solution, eq (76), may be interpreted similarly. The rightward pulse, the first term, prop- 
agates with a decreasing time-dependent effective speed of c(l j—). Due to Heisenberg's inequality, 

4a 

[21] the rightward pulse in F-space gets wider in time. To compensate, so that Parseval's theorem is 
satisfied, the factor in front must decrease more rapidly with time than the corresponding F-space factor. 
The converse observations apply to the leftward pulse. 

We reiterate that it is apparently not possible to go directly from the Fl-space solution, eq (62), to the 
leading-order F-space solution, eq (76), by means of an inverse Laplace transform. The concept of a 
nonunifonnity expansion, as discussed in section 2.1, appears to be meaningful only in F-space and I-space, 
not in F-space or FI-space. Hence we must follow the circuitous route shown in figure 1. 

5. Nonuniform diffusion equation 

5.1. The differential equation 

In the study of the nonuniform diffusion equation, we must at the outset distinguish material diffusion 
from diffusion of a thermal pulse. The equation for material diffusion in one dimension is [23] 

~Jr yiXt t] ~~k- [D{x) ~L~ y{x - i)] = ° {77) 



and the case of interest to us is 



D{x) - D B + D lX (78) 



where D, is small. 

Diffusion of a thermal pulse in one dimension, [24] which corresponds more closely to the Rayleigh- 
Brillomn problem, is characterized by a thermal conductivity X(x) and a heat capacity per unit length o(x). 
The differential equation is [24] 



and the cases of interest are 



X — X„ + Xix 

o — a a + a x x 



(80) 



where X, and a, are small. 

To linear order in nonuniformity, we may represent both equations by the convenient form 

(a - bx) ~f- - c -X + y b Jl- = o «m 

dt dx 1 ' dx ( m > 
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where, for material diffusion, — = D 

a 



— = DJD a (82) 

a 

y = - D„ 
and for diffusion of a thermal pulse, a = oj^ 

b = Oeki - Xoff! (83) 

c = X 2 

7 = X„ X./6. 

The equation thus is characterized by two nonuniformity parameters 6 and 7, and a single small 
parameter 6. 

5.2. Properties of the solution 

Before embarking on a derivation of the solution to eq (81), it is of interest to examine certain properties 
of the solution which follow directly from the differential equation itself. We multiply both sides of eq (81) 
by (1 + bx/a), discard terms quadratic in 6, and integrate over x from —00 to 00. The result is 

3 ["„ y(x, t) dx = (84) 



dt 

so that the norm \y\ considered previously, eq (24), is a constant of the motion. For material diffusion, eq 
(77), the norm is constant for any possible D(x), a consequence of the conservation of diffusing material. 
For thermal diffusion, eq (79), the above norm is not constant in general, but is so to first order in nonuni- 
formity. 

If we consider pulses with unit norm, the central position is given by 

<*> = \:„xy(x,t)dx. (85) 

Multiplying the differential equation by x and integrating yields 

2cb yb 



dt 



<x> = -^- + -L^-. (86) 



Physically, for positive Xi and negative a t the diffusion coefficient effectively increases with x. Here 6 and 
7 are positive, so the average position of the pulse drifts to the right. In the converse situation, the pulse 
drifts to the left. 

A similar calculation yields the time derivative of the mean-square displacement, 

9 , « c / 276 66c \ .__,. 

ir <*> = 2 t + Hr + -*-) <x> • (87) 

and the mean-square width, 

3 <Ax*> =-%-[<x*> - <x> 2 ]=-^£- + -^-<x> (88) 



dt dt 

so that the increase in width with time is independent of 7 
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5.3. Solution 

The solution of the nonuniform diffusion equation closely parallels that of the nonuniform wave equation. 
The differential equation in -Fl-space is 

i a • (89) 

1 ° y(k,0)+-~y(k,0). 



s dk bs 

Again, this linear differential equation in k is solved by the integrating factor technique, with the bound- 
ary condition A, = oo. The result is 



ai -, , i m , 1 d -i 

U T It, V) J 

(90) 



y\k, s) = - \y U [-jt-yXu +k,0)+ —-£-y(u + k, 0) ] 

OS S QK 



T am ci y i 

X exp [— =— + -^~ u(u 2 + 3uk + 3k 2 ) J— u(u + 2k)\ . 

o 5bs is 

Without loss of generality we consider only positive 6. The convolution theorem is used to transform the 
solution to L-space. We are required to add a small positive imaginary part to c. The result is 

y{x, s) = y a{bsicw) \ Q u du exp [-^- + -^ ^-J 



(91) 



X i!L dv y(v, 0)[1 ] exp[-iuv] 

., r cu 3 i ibs , ™ iu , \ 

4cs* 2c 2s J 

Again we make the substitution u = bw and then discard all 0(6 2 ) terms. We next explicitly perform the 
w-integral by means of the substitution w = z 2 and eq (70). It is necessary to assign a small positive imag- 
inary element to a. The result of these operations is 

y{x, s) = ±-(a/csy« Jl dvfa 0)[1 - %1 - f (x + «)]-»'» 
2 a 2a 

■vA A (92) 

X exp H£-(x - v)] exp ~ {(sa/c) in \x - v\ [1 - j- (x + v)] 11 *} . 

At this point it is easily seen that the uniform limit, b — + , recovers the uniform solution, eq (41) with D 
= c/a. 
To determine the f-space solution, we use the relation 

This equation is applicable for real positive a, which requires that 

x + v < -^ . (94) 

o 

Once again, it is built into the mathematics that a solution does not exist for large displacements where the 
original differential equation becomes unphysical. 
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The P-space solution initially is most conveniently represented in kernel form, 

y(*. = J-". y(v, 0) K (x, v;t) dv (95) 

where, to linear order in b, 

K(x,v;t) = ^.{rDjr i2 [l + ±(x - 3v) + -^-(x - v)] 

h , (96) 

X exp {- (x - vf [1 - A( X + t,)]/4Z) t} 

and Z) = c/c is the coefficient of diffusion in the uniform limit. 

However, there is a severe mathematical difficulty with eqs (95) and (96) in their present form; namely, the 
integral of eq (95) does not exist as written. There is an essential difference between the solutions of the un- 
damped wave equation and the diffusion equation. For the undamped wave equation, if the initial pulse is 
strictly localized around the origin, i.e., y{x, 0) = identically for | a; | > w I2, then, for any large distance 
L, y(L, t) remains strictly zero until some finite time of the order L/c. However, as governed by the diffusion 
equation, a similar localized pulse causes a finite displacement at L after any time interval, however small, 
as an infinite tail is instantaneously created of the form exp{ —x 2 ) at large x. Then, since the kernel of eq (96) 
goes as exp( + v 3 ) for large v, the integral of eq (95) is formally infinite after any two successive finite time 
steps. 

Physically, however, it is not possible for a pulse to diffuse over large distances instantaneously; this 
phenomenon should be interpreted as a mathematical artifact. We believe the proper way to interpret eqs 
(95) and (96) physically is to remove this artifact in the following manner. The kernel of the integrand is a 
distorted Gaussian function integrated over an undistorted measure. A transformation is made so that the 
kernel becomes an undistorted Gaussian integrated over a distorted measure. Specifically, we define the 
new measure 

u* =(x -vf[l -i-(x + v)] (97) 

Za 



and 



if = {v - x) [1 - — {x + v)] + 0(6 2 ) . (98) 

4o 



Then, if terms higher than linear order in b are discarded, eq (95) becomes 

y(x, t)=4r (xDjr* 2 jr. du y(x + u[l + -f (2x + «)], 0) 

ri - _^_ - J^Ll e -» 2 '4°o, 
L 4a 2c J 

The major formal change between eq (95) and eq (99) is the upper limit of infinity on the u integration, 
whereas, from eq (97), as v — +<x>, u becomes, formally, infinite and imaginary. The objective of the 
transformation is to suppress the unphysical behavior for x > a/b where the diffusion is mathematically 
negative. Equation (99) is a well-behaved solution to linear order in nonuniformity which connects smoothly 
to the uniform case. 

To check the solution, we first verify that K(x,v;t) satisfies the boundary condition 

K(x,v$) = 5(x - v) (100) 

which is equivalent to 

K(x,v;0) = [x±v] (101) 

and 

£. K(x,v$) dv = 1. (102) 
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Equation (101) is verified by inspection, and eq (102) is seen to be valid in the sense of the transformation 
from eq (96) to eq (99). Finally, it is straightforward to show that the kernel should and does satisfy the dif- 
ferential equation 

c -^- K(x,i;;0 - (a - bx) i- K(x,v;t) - yb£- K(x,v;t) = (103) 

ox* at ox 

through linear order in 6. 

The F-space solution may be derived by expanding y in the integrand of eq (99) through linear order in b, 
rearranging the integrand into sums of terms of the form f(u)g(x — u), and applying the convolution theo- 
rem. The result is 

y(k, t) = y(k, 0) e- D »* 2 ' [1 - i {—■ + -&-) DM + i A Dlk 3 ?] 

A h < 104) 

-[-^f^O^e-^ij-D^t 

This solution satisfies the initial condition by inspection. It is also straightforward to check that it satifies 
the differential equation, eq (81) transformed to F-space, through linear order in 6. 

We proceed to calculate the dynamic structure factor, the time Fourier transform of the intermediate scat- 
tering function F\k,t) (cf. eqs (2-3)), The solution of the initial value problem in principle yields F(k,t) only 
for positive times. We employ the usual convention for a time correlation function that for negative times 
[25] the function is defined to be 

Fflc,t) = F*(k,-t) (105) 

where the asterisk denotes a complex conjugate, in which case 

SM = 2 Re J m e-** F[k,t) dt (106) 

Applying eq (106) to our solution yields the result 

SM> = <H-k, 0)y(k, 0)> { "ff - (JL + J*.) .^J* 

K or + (Dak 2 ) a c kj 2 + (F 



[oj 2 + (D k 2 ff 

+ z a UJi [„. + (zu 2 ) 2 ] 3 / <y( * ,0) dk y{k ' 0)> a [c 2 + (Dry] 1 



(107) 



5.4. Interpretation 

The most important physical features of the solution to the nonuniform diffusion equation may be given a 
simple interpretation as a "Doppler shift" similar to that of the uniform damped wave equation. In par- 
ticular, eq (86) shows that the diffusion pulse, on the average, drifts in the direction of increasing diffusion 
with a velocity c D given by 

2co 70 ..... 

c ° = -jT + "^- (108) 

The uniform damped wave pulse undergoes a similar motion with the wave speed c and an "effective dif- 
fusion rate" of V. By a naive analogy we would expect that the F-space nonuniform solution is 



y(k, t) = y(k, 0) e- D ° k2t e-"° k < 



(109) 
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which, through first order in c D , yields the first two terms of eq (104). By similar reasoning we would obtain 
the dynamic structure factor 

SM = < -n-k, 0)y(k, 0» ( . + ^ m , f (HO) 

which similarly accounts for the first two terms of eq (107) through (c fl ). 

The third term of eq (107) gives a rather complicated change of shape of the "Rayleigh line" not pre- 
dicted by the above heuristic arguments. The term oscillates several times and is qualitatively akin to the 
third frequency derivative of the primary Lorentzian function, so it alters, in leading order, neither the 
mean frequency nor the width. 

The ensemble average of the fourth term in eq (107) vanishes for a uniform system because of invariance 
under a parity transformation. This term is somewhat difficult to interpret, but, when averaged over an 
appropriate ensemble, is presumably no more than first order in nonuniformity. Hence the entire term pre- 
sumably yields a contribution second order in nonuniformity and may be neglected. 



6. Nonuniform damped wave equation 

6. 1 . The differential equation 

We first seek the appropriate one-dimensional nonuniform analog to the equation, first derived by Stokes, 
for the propagation of a plane sound wave in an isotropic viscous medium, [26] 

^-"-^--T^-^"^- ' (111) 

cf. eq (46), where q is the density, f the displacement, x the effective bulk modulus, and /* the viscosity. This 
corresponds to the neglect of coupling between density and temperature fluctuations in the Rayleigh- 
Brillouin scattering problem. [11] 

The following model, though somewhat contrived, appears to be the one-dimensional analog to the above 
with the appropriate properties. The viscous term corresponds in three dimensions to an off-diagonal term 
in the pressure tensor, which, crudely speaking, represents a torque per unit volume. In one dimension, this 
would be a torque per unit length. 

We consider a stretched string with position-dependent mass per unit length q(x). We imagine that, as the 
string vibrates transversely, it traverses a viscous medium which exerts a torque AL on a length element A*. 
On each element Ax there is a force F[x), perpendicular to the string, acting at one end and an equal and 
opposite force acting on the other end. Thus 

AL = [F[x) + F[x + Ax)] Ax . (112) 

We further assume that the force is proportional to the spatial derivative of the velocity of the segment, 
the analog of the three-dimensional velocity gradient 

JW-^W^L-fL. (.13) 

where fi{x) is analogous to viscosity and may be position-dependent. 

Next we generalize the standard derivation [27] of the wave equation for a stretched string under tension 
Thy the inclusion of this new force. As usual it is assumed that the angle between a vibrating segment of the 
string and the equilibrium direction is sufficiently small that its sine and tangent are equivalent. Applying 
Newton's second law to the segment of the string from x x to xi yields 
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or, since x l and x 2 are arbitrary, 

This is the particular nonuniform damped wave equation to be analyzed here. We are interested in the 
special case 

q(x) = a + bx (116) 

liipc) = ft + 0x (117) 

where b and are small. 

6.2. Solution 

The solution of the nonuniform damped wave equation closely parallels those of the undamped wave and 
diffusion equations. The equation in FL-space is 

as z y + ibs 2 -|f- + TBy + IpsWy + 20iksy + 2/3jPs 4£- = as y(k, 0) 



+ ibs-^-y(k, 0) + 2/i* 2 y(fc, 0) + 20ik flfc 0) + 2#it 2 -^- y(it, 0) 



As before, this linear differential equation may be solved by the integrating factor technique. [20] 
However, at this point we anticipate a mathematical difficulty. In previous examples, where the linear dif- 
ferential equation had the form 

-f*- +./W r = «(*> (U9) 

the variable y(&) had the form of a quadratic polynomial in k. This property made possible the subsequent 
inverse Fourier transformation (from FL-space to L-space) by means of eq (65). Upon examination of eq 
(118), we note that J[k) is quadratic if and only if = 0. We shall therefore set = in the remainder of this 
section. This corresponds to uniform "viscosity" with nonuniform density, or a particular correlated case of 
nonuniform wave speed and wave damping. 



Again we use the boundary condition k t = oo, and make the substitution k' = u + k [cf. eqs (61)-(62)]. 
The FL-space solution is 

y{k, s) = I> [-fj-r(" + k, 0) + 2 'Hfr+ff y {u + k , 0) 

- \l k T\u + k, 0)] exp[- -SSL -il±^L u{ui + Uu + 3^)] (1 
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Without loss of generality, we consider only positive b. The convolution theorem is then used to derive the 
L-space solution. Once again, it is necessary to assign a small negative imaginary component to T. The 
result is 



yix, s) = f- [irbi {T + 2ns)]- U2 \'du u xn exp [ - 



lau 



2 ^.^ ^ n Jo ~~~ ~. rL b 



(r+2^)m 3 i r- r bv . 2 M a 2 . „i 

3b?—* 1- dv W v > W 1 + — > - -^-W^' 0) J 

X «*,(-»,,) e>;> ( T [ ^ 2u(x -v)+ {T+2fis)u J ) 



(121) 



As in the undamped wave solution, we make the substition u = bw and discard terms explicitly 0(b 2 ). We 
then let w = z 2 and apply eq (70), after assigning a small positive imaginary part to a. The result is 

SM = Via* (T+2psF* \% dv [yiv, 0)[1 + -^-] - -^— ^-^ 0)] 

[l +-^~(x + v)[- XA exp {-,|x-»|[l + -t-( x + v )]»[ a l(T + 2ps)]»} . (122) 

In the uniform limit, b — + , we recover eq (52), with r = /i/a and c = (T/a) x/i . 

At this point we encounter the problem of dispersion discussed in section 2.3. Since there is no general 
P-space solution to the uniform damped wave equation in terms of elementary functions, we cannot expect 
to find an explicit P-space solution to the nonuniform equation either (to all orders in dispersion). A formal 
expression may be written, 

*M = V4 L°> [1 + -^-(* + vV* [y(v, 0)(1 + ~^-)Xr l {(c 2 + 2TsF* 

exp [s\x - v\ [1 + -L-( x + v )]*(c* + 21^)"*]} - ZT-^yiv, 0)£~ l {s~ l (123) 

(c 2 + 2r 5 )-'^ exp[s\x. - v\ [1 + -4—(x + v)P (c 2 + 2^)"'^]}] 

where-£ -1 is the inverse Laplace transform operator. However, this form is not particularly illuminating 
physically. 

To display a more transparent solution, we assume that there exists an operator^" 1 which simultaneously 
performs an inverse Laplace transformation and projects out the leading term in a dispersion expansion. 
For example, on comparing eqs (52) and (55), we see that such an operator should effect the transformation 

I' 1 {(c 2 + 2rs)-» exp[-s\x - v\ (c 2 + ZTsY*]} 

= 1 /2(7rr*)"Me~ (x_r+c ' ), ' 4r ' + e-b-'-t'W]. (124) 

We further assume that-^ 1 " 1 obeys the change of scale property of Laplace transforms, namely, if 

£ x {As)} = F(t) (125) 

then 

-f ■ {*»>} = >(^) . (126) 

at least providing a depends only on position and nonuniformity parameters. The above is, in essence, a 
statement of the independence of nonuniformity and dispersion expansions. 
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The second term of eq (116) is of higher order in dispersion than the first term and may be neglected. On 
replacing-^" 1 by-?" 1 in eq (123) and using eqs (124)— (126), we obtain the solution 

Ax, t) = K(*nr* O y(v, 0)[1 +^][1 +^x + «)]-» 

[exp { - [ct - (1 + -^- [x + v]p (x - vfflWt} (127) 

+ exp { -[ct + (1 + -^- [x + v])» (x - v)] 2 I4Tt}] 

which bears a strong similarity to the nonuniform diffusion solution, eqs (95) — (96), for 7 = 0. 

As in the case of diffusion, we encounter problems with divergences in the integrals due to mathematical* 
ly instantaneous propagation of the pulse to points at large distances from the origin, where the original dif- 
ferential equation makes no physical sense. We again surmount the difficulty by transforming to a distorted 
measure in which the Gaussians of the kernel are undistorted. The new measure is 

u=(x-v)[l +-^-(x + v)}» (128) 

Then a calculation, discarding terms higher than linear order in b, analogous to that of diffusion yields 
y{x, t) = Vtbrrt)-* II du y[x ~ u(l - -^- [2x - u]), 0] 

(1 — ^-) [e-fa-o-ur, + e -(.. + c0'/4r,]. (129 ) 

4o 

Since 

Urn ti (icTt)-* e -(.±o-/4r« = % u ±ct) (130) 

it is easily checked that, in the limit of zero damping, the nonuniform undamped wave equation solution, eq 
(75), is recovered. 

We recall that in the present problem there are two small expansion parameters, namely, nonuniformity 
and dispersion. Since the latter is defined in terms of .F-space variables, there is no natural means of check- 
ing that the P-space solution satisfies the differential equation. 

The f-space solution may be derived by expanding y in the integrand of eq (129) through linear order in 
b, rearranging the integrand into sums of terms of the iormJ[u + ct)g(x — u), and repeatedly applying the 
convolution theorem. The end result is rather lengthy and is best expressed as the sum of a rightward and a 
leftward pulse. 

y(kt)=y r (k,t) + y,(k,t) (131) 

where 



frikt) = y{k, 0) e-™ 2 ' e- Uk ' 
+ [-JL- yik 0)] e-" 1 ' *-■* [-A- Ml - -A- cht] 



(132) 



and y t is obtained from y, by replacing c by — c throughout. Since we only desire the solution through 
linear order in dispersion, the term above containing T 2 may be discarded. 

450 



The F-space transformation of the original differential equation, eqs (115-117) with /3 = 0, is 

■£ + T i -n--g- + *vf + 2 « , -£ = » (133) 

At this point it can be checked that the solution, eq (132), satisfies the differential equation, eq (133), if terms 
0(6 2 ), second order in nonuniformity, and 0(r 2 ), second order in dispersion, are neglected. 

The solution satisfies the initial value condition by inspection, but does not satisfy the condition on the 
initial time derivative. The latter, however, is not a failing of the nonuniformity analysis but, rather, a 
feature of the truncation of the dispersion expansion of the uniform damped wave solution. In particular, eq 
(51) is the correct solution to all orders in dispersion, whereas eq (54) satisfies the differential equation and 
initial value condition through first order in dispersion, but not the initial time derivative condition. A fully 
consistent first-order truncation of the dispersion expansion which includes the time derivative condition 
would retain the second term of eq (51) and replace (c 2 — r 2 £ 2 ) I/2 by c throughout. 

A similar difficulty appears in the standard treatment of Rayleigh-Brillouin scattering; compare, for 
example, eq (27) of ref. [1 1] with eq (97) of ref. [28]. Although the extra term required to make the initial time 
derivative vanish is in practice often neglected, for the sake of a fully consistent solution through first order 
in both dispersion and nonuniformity we shall find and include such terms here. 

We note that eq (131) remains a solution through the required order if terms are added which are propor- 
tional to either V times the solution with r — in the brackets of eq (132), or b times the solution with b — • 
in the same brackets. It is straightforward to set up a trial solution with variable coefficients for the re- 
quired additional terms, and solve for the coefficients subject to the constraints that the initial value vanish 
and the initial time derivative be equal and opposite to that of eqs (131)— (132). To obtain the desired result, 
a term 



AyXkt) = y{k, 0) e- Tk2t e 



Tk2t.-ickt 



L c v 2 8a 8a ' 8ca J 



(134) 
+ [-4-y[k, 0)] e -r*V'* [_Tkb_ _ ilM] 
ok 4ca 4a 

must be added to y r (k,t\ and a term Ay ( (k, t), eq (134) with c replaced by -c throughout, must be added to 
y ( (k, t). We could, in a similar manner, impose a nonzero initial time derivative proportional to b and deter- 
mined from external considerations, but will not do so here. 

In the construction of the dynamic structure factor we will, for reasons decribed in section 5.4, neglect 

terms involving ( „ y(k,0).) The resulting dynamic structure factor may be written as 

S(k,u>) = SAk,u) + Sfcu) (135) 



where 



SXM = <y(-k,0)y(k,0)> X [(1 + -|P-) m 



Sea ' (u + ckf +(rk 2 f 

5b p 2 p a + ck 

c (u + ckf + (fit 2 ) 2 ' 2a [(oi + ckf + (r* 2 ) 2 ] 2 



_Tk cu + ck + 5b p 2 p Q3 + ck 



(136) 
_ Sbc_ (Tk 2 ? - (u + ckf Zbc_ {V „„ {Tk 2 ? -3(a) + ckf 

4a [(u + ckf + (rit 2 ) 2 ] 2 2a K F [(u + ckf + (r* 2 ) 2 ] 3 

b<?k 3(u + ck)(rVY - {u + ckf i 
2a [(oi + ckf + (rP) 2 ] 3 



and 5/A,fa)) is given by eq (136) with c replaced by -c throughout. 
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6.3. Interpretation 

In order to develop a physical understanding of the nonuniform damped wave equation solution, we first 
note that the effects of nonuniform speed and nonuniform damping are not independent and, in the special 
case considered here, tend to cancel each other. According to eqs. (115-117) with b > and |3 = 0, both the 
wave speed and attenuation decrease to the right and increase to the left. Hence we expect the rightward 
pulse to undergo less damping, effectively, than in the uniform case. However, as has been seen explicitly 
with the nonuniform undamped wave solution, the decreasing wave speed attenuates the rightward pulse in 
a different sense by lowering its time-dependent integrated norm [eq (24)]. 

We have been careful to retain all terms linear in the nonuniformity and/or dispersion parameters. 
However, some terms in the end give much smaller contributions to S(k,u>) than others. For example, the 
conditions from table I on the validity of 5(£,a>) together with the smallness of the dispersion parameter 
imply the inequalities 

< < -^- < < 1 . (137) 



ah 



The rightward pulse (with our sign conventions) corresponds to the Brillouin line (Stokes line) centered at 
oj = — ch. The additional factor STbl&ca multiplying the primary Brillouin line [first term in eq (136)] is by 
(137) at least three orders of magnitude less than unity, and thus utterly negligible. 

The remaining terms in eq (136) are conveniently classified quantitatively, to leading order, in terms of 
the shift in mean frequency or change in width they cause in the Brillouin line. Since the more commonly 
used definitions of mean position or width of a peak involve integrals which diverge for the Lorentzian func- 
tion, we employ the following alternate definitions, for mean frequency, 

w = II « [/Tco)] 2 dul J_-„ [/fa)] 2 do, (138) 

and for width of the line, 

(Au) 2 = II (co - Hf \JW)Y do,l II [/Ico)] 2 da, . (139) 

For the primary Brillouin line eq (139) gives, as usually defined, Aco = Th 1 . The precise definitions are not 
important but the above ones serve as convenient measures of changes in the spectrum for purposes of 
discussion. 

The second term of eq (136) is antisymmetric about a, = — ch and corresponds to the well-known first- 
order dispersion term in the Rayleigh-Brillouin spectrum [28] in the limit that the specific heat ratio 
approaches unity. To leading order in dispersion this term causes a shift of {2T 2 h 3 lc) to the right in the mean 
frequency of the line. 

The third term of eq (136) is also antisymmetric about oj = — ch. It has the effect of shifting the mean fre- 
quency to the right by (5bThla). This represents a very small shift, which by (137) is less than that of the 
previous dispersion term which is itself usually neglected in practical applications. 

Upon comparison with the undamped wave solution, the fourth term of eq (136) is seen to originate from 
the decreasing wave speed or increasing inertia in the rightward direction. This term is symmetric about co 
= — ch and, to leading order, increases the width of the line by Zbc/4a. An increase in width corresponds 
physically to a larger effective attenuation coefficient. 

The fifth term of eq (136) is symmetric about oj = — ch and, by the definition (139), results in a decrease in 
width of 36c/4a. Thus the fourth and fifth terms effectively, if not identically, cancel each other. One possi- 
ble physical interpretation of the latter term is as follows. The Fourier component of wavenumber k of the 
rightward pulse has a lifetime of (Th 1 )' 1 , and therefore a mean path length cITh 2 . Equations (115-117) imply 
a distance-dependent attenuation factor 



T(x) = 1X1 - ^) (140) 
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from which we infer that the effective attenuation coefficient of the kth Fourier component of the rightward 
pulse is 

r. // = r-/ J -^- (i4i) 

where/! is a positive number of order unity, and the resulting width of the Brillouin line is r, /f k 2 . The full 
analysis then gives /i = 3 A- 

Hence the near-cancellation of the previous two terms may be an artifact of setting /3 = in eq (117). For 
a more general solution where, for example, wave speed and damping increase in opposite directions, the 
corresponding new terms might reinforce each other. In any case, measurement of such a width change 
would be difficult as it would require a fractional precision in the width measurement equal to the fractional 
change in attenuation over the path length of the monitored Fourier component. 

Finally, the last term in eq (136) is highly oscillatory and qualitatively similar to the third derivative of the 
Lorentzian. Like its counterpart in the diffusion solution, it would be very difficult to detect and may be 
neglected. 

Similar observations, with appropriate changes of sign in some cases, apply to S,(i,aj), the anti-Stokes line 
centered at oj = ck resulting from the leftward pulse. 



7. Summary 



As an initial step toward the understanding of Rayleigh-Brillouin scattering in a temperature gradient, 
which involves nonuniform, three-dimensional coupled damped wave and diffusion equations, we have 
analyzed in detail the initial value problem for a set of prototype one-dimensional differential equations 
with linearly nonuniform coefficients, to first order in nonuniformity. 

A question which first arose was whether explicit boundaries were necessary for the system, since at 
large distances the differential equations become unphysical. For the models studied, it is possible to derive 
meaningful solutions of the initial value problem without such boundaries for sufficiently narrow initial 
pulses and sufficiently short times. Also, if the equation has a damping mechanism, it is possible to construct 
a meaningful dynamic structure factor for values of k whose Fourier components are damped sufficiently 
rapidly. Conditions for the validity of solutions are listed in table 1. 

Our procedure has been first to Fourier transform in space and Laplace transform in time. The resulting 
equations, algebraic in the uniform case, become first-order differential equations in k for the linear non- 
uniform case. An expansion through quadratic order in nonuniformity would lead to a second-order differ- 
ential equation in k, for which there is no general method of solution. 

We have seen that the condition of small nonuniformity cannot be represented by any dimensionless 
parameter in F-space or FZ^space. Therefore, to construct the time evolution of the pulse and, ultimately, 
the dynamic structure factor, we must follow the circuitous path indicated in figure 1. The ability to perform 
the initial inverse Fourier transform restricts the range of solvability of the linear differential equation in k. 
One of our initial goals, the solution of the nonuniform damped wave equation with arbitrary, independent 
variation of the nonuniformity in both wave speed and damping, for this reason has not been realized. 

We have derived explicit solutions to the initial value problem for the nonuniform undamped wave, diffu- 
sion, and (for a special case) damped wave equations. These solutions display properties expected on intui- 
tive grounds, and reduce properly to their respective uniform solutions. 

Throughout the derivations, certain mathematical liberties have been taken. We have assigned small 
imaginary parts to intrinsically real coefficients, although these imaginary components go to zero at the end 
of the calculation. For the diffusion and damped wave equations, we have altered the measure of integration 
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in the final solution to suppress singularities at large distance. In the damped wave equation, we have intro- 
duced and postulated properties of an operator which simultaneously performs an inverse Laplace trans- 
formation and projects out the leading term in a dispersion expansion. Attempts to justify these liberties 
have been based on pragmatism and physical intuition rather than mathematical rigor, although we conjec- 
ture that such steps could (probably) be justified mathematically. Of course, an a posteriori justification is 
that our solutions explicitly satisfy the differential equations through linear order in nonuniformity, and 
that our initial conditions specify a unique solution. 

Since the present formalism neglects coupling between wave and diffusion modes, three-dimensional 
effects, finite initial time derivatives and position-dependent rates of generation of fluctuations, one should 
be cautious in drawing conclusions from these results about experimental Rayleigh-Brillouin scattering 
from a fluid in a temperature gradient. However, the results suggest that further studies of possible changes 
in the Rayleigh line could be fruitful. Kirkpatrick, et al. [6] have noted that certain modifications of the 
Brillouin spectrum vanish when an average over a finite volume of the fluid is taken. But according to the 
simple and intuitive arguments of section 5.2, all heat modes within a given volume should experience a 
drift velocity in the direction of increasing thermal diffusivity. Our one-dimensional results suggest that the 
Rayleigh line is thus Doppler shifted. Detection of such a shift, if it persists in the coupled three- 
dimensional problem, depends on precision calibration of the incoming laser frequency and precision 
measurement of the mean frequency of the Rayleigh line. The shift would be most pronounced for systems 
whose thermal conductivity changes rapidly with temperature and for large values of k. By contrast, the 
asymmetry in the size of the Brillouin peaks predicted by previous authors, [1-10] and not considered here, 
is generally agreed to be inversely proportional to k 2 , and experiments to detect such an effect [14] have thus 
been performed at very small values of k. 

The author thanks R. D. Mountain for elucidating many aspects of light scattering theory and for guid- 
ance and encouragement. He thanks H. J. Raveche and the staff of the Thermophysics Division, National 
Bureau of Standards, Washington, DC for their hospitality during the author's visit in 1979, when much of 
this work was completed. He also thanks S. W. Haan, R. J. Rubin and J. R. Dorfman for valuable 
suggestions. 



8. References 

According to conversations with J. R. Dorfman and G. van der Zwan, the authors of refs. [1-10] have now 
formed a consensus that the correct leading-order spectrum is that given in refs. [6, 9, and 10]. 

[1] Procaccia, I.; Ronis, D.; Oppenheim, I. Light scattering from nonequilibrium states: the implication of broken time-reversal sym- 
metry. Phys. Rev. Lett. 42(5): 287-291, 1979 January 29. 
[2] Procaccia, I.; Ronis, D.; Collins, M. A.; Ross, J.; Oppenheim, I. Statistical mechanics of stationary states. I. Formal theory. Phys. 

Rev. A 19(3): 1290-1306; 1979 March. 
[3] Ronis, D.; Procaccia, I.; Oppenheim, I. Statistical mechanics of stationary states. II. Application to low-density systems. Phys. Rev. 

A 19(3): 1307-1323; 1979 March. 
[4] Ronis, D.; Procaccia, I.; Oppenheim, I. Statistical mechanics of stationary states. III. Fluctuations in dense fluids with applications 

to light scattering. Phys. Rev. A 19(3): 1324-1339; 1979 March. 
[5] Kirkpatrick, T.; Cohen, E. G. D.; Dorfman, J. R. Kinetic theory of light scattering from a fluid not in equilibrium. Phys. Rev. Lett. 

42(14): 862-865; 1979 April 2. 
[6] Kirkpatrick, T.; Cohen, E. G. D.; Dorfman, J. R. Hydrodynamic theory of light scattering from a fluid in a nonequilibrium steady 

state. Phys. Rev. Lett. 44(7): 472-475; 1980 February 18. 
[7] Van der Zwan, G.; Mazur, P. Light scattering from a fluid with a stationary temperature gradient. Phys. Lett. 75A(5): 370-371; 

1980 February 4. 
[8] Van der Zwan, G. On some aspects of fluctuating hydrodynamics. Thesis. Ryksuniversiteit te Leiden; 1980. 
[9] Tremblay, A.-M. S.; Siggia, E. D.; Arai, M. R. Fluctuations about hydrodynamic nonequilibrium steady states. Phys. Lett. 76A(1): 

57-60; 1980 March 3. 
[10] Ronis, D.; Putterman, S. Simplified hydrodynamic therory of nonlocal stationary state fluctuations. Phys. Rev. A 22: 773-777; 

1980 August. 
[1 1] Mountain, R. D. Spectral distribution of scattered light in a simple fluid. Rev. Mod. Phys. 38(1): 205-214; 1966 January. 
[12] Komarov, L. I.; Fisher, I. Z. Theory of Rayleigh scattering of light in liquids. Soviet Phys. J.E.T.P. 16(5): 1358-1361; 1963 May. 
[13] Pecora, R. Doppler shifts in light scattering from pure liquids and polymer solutions. J. Chem. Phys. 40(6): 1604-1614; 1964 

March 15. 

454 



[14] Beysens, D.; Garrabos, Y.; Zalczer, G. Experimental evidence for Brillouin asymmetry induced by a temperature gradient. Phys. 

Rev. Lett. 45(6): 403-406; 1980 August 11. 
[15] Brownell, J. H. Wave number space analysis of propagation in nonuniform media. Am. J. Phys. 41(2): 207-212; 1973 February. 
[16] Clark, N- A. Inelastic light scattering from density fluctuations in dilute gases. The kinetic-hydrodynamic transition in a 

monatomic gas. Phys. Rev. A 12(1): 232-244; 1975 July. 
[17] Landau, L.; Placzek, G. Structure of the undisplaced scattering line. Physik Z. Sowjetunion 5: 172; 1934. 
[18] Reif, F. Statistical and thermal physics. New York: McGraw-Hill; 1965; p. 301. 
[19] Mountain, R. D. Generalized hydrodynamics. Adv. Mol. Relax. Proc. 9: 225-291; 1976. 
[20] Mathews, J.; Walker, R. L. Mathematical methods of physics. New York: Benjamin; 1965; p. 3. 
[21] Dym, H.; McKean, H. P. Fourier series and integrals. New York: Academic Press; 1972; p. 116. 
[22] Gradshteyn, I. S.; Ryzhik, I. M. Tables of integrals, series and products. New York: Academic Press; 1965; p. 307. 
[23] Richardson, L. F. Atmospheric diffusion shown on a distance-neighbour graph, Proc. Roy. Soc. London Series A 110(756): 

709-737; 1926 April. 
[24] Berg, P. W.; McGregor, J. L. chapter 3 in Elementary partial differential equations. San Francisco: Holden-Day; 1966; 29-33. 
[25] McQuarrie, D. A. Statistical mechanics. New York: Harper and Row; 1976; p. 563. 
[26] Crandall, I. B. Theory of vibrating systems and sound. New York: Van Nostrand; 1926; p. 110. 

[27] Berg, P. W.; McGregor, J. L. chapter 8 in Elementary partial differential equations. San Francisco: Holden-Day; 1966; p. 181-187. 
[28] Mountain, R. D. Liquids: dynamics of liquid structure. Chem. Rubber Co. Crit. Revs. Solid State 1(1): 5-46; 1970. 



455 



